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A lower bound is derived for the boundary entropy s = In g of a 1+ld quantum critical system 
with boundary under the conditions c > 1 on the bulk conformal central charge and Ai > (c— 1)/12 



ON on the most relevant bulk scaling dimension. This is the first general restriction on the possible 



values of g for bulk critical systems with c > 1. 
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A 1+ld quantum critical system is described by a 2d conformal field theory (the bulk 
CFT). A critical boundary is described by a conformally invariant boundary condition on the 
bulk CFT. The combination — a bulk CFT with a conformally invariant boundary condition 
- is a boundary CFTflj. Critical junctions in critical quantum circuits are described by 
boundary CFT's. For an iV-wire junction, the bulk CFT is the iV-fold product of the 
CFT's describing the individual wires. The critical junction is described by a conformal 
boundary condition for the product CFT. In string theory, branes in spacetime are described 
by conformal boundary conditions on the string worldsheet. 

Affleck and Ludwig[2] defined a number g for each boundary CFT — the universal non- 
integer ground state degeneracy. The entropy localized in the boundary — the boundary 
entropy — is In g. It is defined as the total entropy of the system minus the bulk entropy 
7rcL/6/3, which is proportional to the length L in the limit of large L. The coefficient of 
L is determined by conformal invariance, in terms of the inverse temperature /3 and the 
conformal central charge c of the bulk CFT. 

For c < 1, there is a complete classification of all possible boundary CFT's[3j. There 
is also a complete classification of conformal boundary conditions for the c = 1 gaussian 
model[lH6]. Until now, no limitations have been known on the possible values that g can 
take for any other c > 1 bulk systems. 

For non-critical boundary conditions in a bulk CFT, the boundary entropy s is defined 
in the same way by subtracting the universal bulk entropy from the total entropy. Now s 
depends on the temperature. Under a change of the thermal length scale j3 the effective 
boundary condition evolves along the boundary renormalization group flow (the boundary 
RG flow). The bulk system, being scale invariant, stays the same. A fixed point of the 
boundary RG flow is a boundary CFT. At a fixed point s = Ing. It is not obvious that s 
decreases with decreasing temperature — that the second law of thermodynamics applies 
to the boundary — because of the subtraction of bulk entropy in the definition of s. In 
fact, the boundary entropy s does decrease along the boundary RG flow, so it decreases 
with temperature [7]. The result is actually stronger: the boundary RG beta function is the 
gradient of the function s on the space of boundary conditions. All that is missing to control 
the asymptotic low temperature behavior is a lower bound on s. Such a lower bound would 
be an analogue of the third law of thermodynamics. Again, the existence of a lower bound 
on s is non-obvious because of the subtraction of the bulk entropy. Unsuccessful attempts 
have been made to prove that s is bounded below [8.]. Without a lower bound, we cannot 
exclude the possibility that s might decrease to — oo as the temperature drops to zero. 

Here we prove a lower bound g > gs(c, Ai) that applies to any c > 1 bulk system that 
has Ai > (c — 1)/12 where Ai is the most relevant bulk scaling dimension. The proof 
assumes nothing about the boundary condition besides criticality and unitarity. The bound 
does not imply a boundary third law of thermodynamics, since it applies only to critical 
boundary conditions. It does imply that a non-critical boundary with entropy s below the 
bound cannot flow to a critical boundary condition at zero temperature. If such a system 
exists, its boundary entropy must necessarily decrease without limit towards s = — oo at 
zero temperature. 

One of us has argued that critical quantum circuits are natural physical systems for 
asymptotically large scale quantum computers [9]. The quantum wires should be critical in 
the bulk, so that the low-energy excitations are protected against microscopic fluctuations 
by universality (the RG), and travel at uniform speed. The processing elements are to be 
the circuit junctions. A junction can be considered as a boundary condition on the CFT 
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describing the independent wires entering it. A lower bound on \ng leads to an upper 
bound on the information capacity of the junction, giving a general constraint on the design 
of critical quantum circuits. 

In string theory, g is the brane tension. The lower bound on the brane tension might 
be useful once it is extended to superconformal boundary CFT's and if the condition A x > 
(c — 1)/12 can be relaxed. 

The modular duality formula for a boundary CFT is [TO] 



tr exp (-/3H bdry ) = (B\ exp(-2irH bulk / (3) \B). 

On the left is the thermodynamic partition function Zl{(3) at inverse temperature (3 for a 
finite segment of the system of length L — 1. The boundary conditions at the two ends of the 
segment are the same. The hamiltonian is H b dry The Hilbert space is called the boundary 
sector (in string theory, the open string sector). In the euclidean space-time interpretation, 
Zl{(3) is the partition function of a finite 2-d cylinder with length L and euclidean time 
periodic with period (3. The right hand side is obtained by re-interpreting L as euclidean 
time and (3 as the length of a circle or, by scale invariance, euclidean time 2ttL/(3 and a 
spatial circle of length 2tt. The hamiltonian for the circle is H bu ik- The Hilbert space of the 
bulk system on the circle is called the bulk sector (the closed string sector). The boundary 
condition on each end of the cylinder is described by a bulk state \B). The modular duality 
formula states that the partition function depends only on the 2-d geometry, so the two 
quantum mechanical interpretations give the same result. 

Conformal invariance implies that each side of the duality formula can be expressed as a 
sum over the characters of the irreducible unitary representations of the Virasoro algebra. 
For c > 1 (we consider the case c = 1 separately below) the duality formula becomes 



/?) + E X* ) = S 2 Xo(V/3) + E b " XA fc /a(«/0) 



where the characters Xh(0) are given by 



f (ff\-I<T , 0--q), h = 



q = e~^, r^W^na-O, 7 

n=l 



24 



The character XofyP) is the contribution to the partition function from the boundary sector 
representation that contains the ground state, whose energy is — 27rc/24. The characters 
Xhj{if3) are the contributions from the representations with lowest energies 2ir(hj — c/24). 
Unitarity and uniqueness of the ground state imply all hj > 0. The boundary scaling fields 
are in one-to-one correspondence with the energy eigenstates in the boundary sector, via 
radial quantization. A primary boundary field of scaling dimension hj corresponds to the 
lowest energy state in the representation labelled by j. The bulk scaling fields are in one- 
to-one correspondence with the energy eigenstates in the bulk sector. The terms on the 
right side of the duality formula come from the closed sector representations whose lowest 
energy states correspond to the spin-zero primary scaling fields whose scaling dimensions 
are < Aj < A2 < • • • . The numbers g and b k characterize and completely determine the 
conformally invariant boundary state \B). 
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Rattazzi, Rychkov, Tonni, and Vichi[TT] developed the linear functional method for de- 
riving bounds on the low-lying scaling dimensions of conformal field theories from crossing 
formulas for correlation functions. Hellerman[T2] showed that the same method could be 
applied to the modular duality formula for the bulk partition function of a 2d CFT to obtain 
an upper bound on the dimension of the lowest non-trivial primary field, and, with one of 
us, to obtain bounds on state degeneracies[13j. Here we apply the linear functional method 
to the modular duality formula for boundary CFT to derive a lower bound on g. 

We want a bound on g that depends only on properties of the bulk system so it will 
apply to all possible critical boundary conditions for a given bulk critical system. The 
derivation should use only universal facts about the boundary condition: the uniqueness 
of the boundary sector ground state and the positivity of the scaling dimensions hj, which 
follows from unitarity. Otherwise, nothing should be assumed about the numbers hj or b^. 

We start by multiplying both sides of the duality formula by r)(i/3) = (3~ l / 2 ri{i/ (3) to get 

/O + = //O + X>*/A fc (1) 
3 k 

where 

/a - I /3 -i/2~- 7+ a/2 ) A>0; Q-e 

Then we apply a linear functional — a distribution p((3) — to both sides of equation ([!]), 
giving 

(p, / ) + h 3 ) = g 2 (p, /o) + J2 b t (P. /a J 

3 k 

where (p, F) = J °° df3 p((3) F((3) . If we can choose p(/3) so that 

(p,A)>0, Vh>0 (2) 
(p,/a)<0, VA>A x (3) 



then we get an inequality 
Next, using the identity 



g 2 (P, fa) > (P, fo) ■ (4) 



-l/2~-y+A/2 = / dye -*^Wa=*y (5) 



we see that condition (|2j) implies (p, fo) > so we have a lower bound on g, 

2 \ 2 r 1 \Pi f°) lr\ 

g >gB[p\ = ~, jt • (6) 
{p, fo) 

Maximizing over all distributions p(/3) satisfying conditions ^ and (|3]), we obtain the opti- 
mal bound 

g 2 > g 2 B (c, Ai) = max p ^|[p] . (7) 

It is not obvious that there exists any distribution p(/3) satisfying both conditions ^ and 
Using identity ([5]), condition ^ requires 

/ dy (p, f 1+y 2 /2 ) cos fay ^A, - 2jj < . (8) 
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If Ai < 27 this is incompatible with condition (J5]). So the linear functional method can give 
a bound only if A 1 > (c - 1)/12. 

The next step is to approximate the space of distributions by distributions of the form 
(p, F) = VF(f3) where V is an N th order differential operator in (3. A bound g%(c, A 1; N, (3) 
is obtained by taking the maximum in equation ([7]) over the differential operators of order 
N. The bound can only improve as N increases. The partition function is real analytic in (3 
so we can expect the limit N — > 00 to exhaust the space of linear functionals for any choice 
of f3, giving the optimal bound g%(c, Ai) = liniA^oo g%(c, Ai, N, (3). We stop here at N — 1, 
contenting ourselves with finding any bound at all. Elsewhere we will use the numerical 
techniques of [H] (semi-definite programming) to approximate the optimal bound from the 
linear functional method. 

For N = 1, we write the general first order operator 

fid 

V = a + ai - — 



2tt<9/3 

Condition g is do, ai > 0. There is no bound if a\ = 0, and the bound does not change if 
we scale V, so we might as well set a% — 1. Condition ^ then becomes 

These conditions require A\((3) > which cannot be satisfied for any value of (3 if Ai — 27 < 
0, so to get a bound we have to assume A x > 27, the necessity of which we have already 
seen from the general analysis. The bound ^ is 

where 

AM^^q-^l^, A 3 (/3) 



1 — q 1 — q' 



4vr/3 f3 2 /3 2 l-q 

Since ^2,3,4^) > 0, the highest bound is obtained when ao takes its maximum value Ai({3), 

so 

^■ A ^ = A ^W^$ry (10) 

The bound is empty unless A\{f3) — A 3 ((3) > 0, which is stronger than Ai((3) > 0, so 

A 1 (^)-AM>0 (11) 

is the only condition we need to impose to get a bound. 

At this point there is no reason to stick to one particular value of (3. The dependence on 
f3 will disappear as iV — > 00 but for finite iV we can sample a larger subspace of distributions 
if we vary (3. The best bound that can be obtained with a first-order T> is 

0i(c,Ai,l) = max0i(c,Ai,l,/9) (12) 
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where the maximum is taken over all /3 satisfying condition (11). There is a unique positive 



solution f3\ of Ai(f3i) — A$(f}\) = and condition (11) is equivalent to < /3 < So for 
Aj > 27 there is a lower bound 

g 2 >g 2 B (c, A 1 ,l) 

with 

'■^■■"-a^ Htw ' (13) 

There is no analytic expression for the N = 1 bound, but it can be calculated numerically 
for any given value of c and Ai. In general, the detailed form of the N = 1 bound as a 
function of c and Ai is not particularly interesting since it is not even the optimal linear 
functional bound. At this stage, we are only interested in the fact that there is any lower 
bound on g. 




FIG. 1. The N = 1 bound for c = 1 compared to the minimum value of g 2 for the c = 1 gaussian 
model. The comparison is extended past the maximal value Ai = 1/2 by interpreting Ai as the 
lowest dimension of a primary occurring in the boundary state. 



For c = 1, there are degenerate Virasoro representations that do not occur for c > 1. 
Equation ([!]) holds with the modification that, for integers n > 1, f n 2 = g~ 7+n2 (l — q 2n+1 ), 
f2n 2 — (3~ l l 2 q~ 1+n (1 — q 2n+l ). As before, we apply a first order differential operator with 
appropriate positivity conditions to get a lower bound on g that depends on /3, then we 
maximize over 0. We omit the calculations. The result is shown in Figure [TJ Included for 
comparison is the smallest value of g 2 for the c = 1 gaussian model. The N = 1 bound is 
moderately good except when Ax ~ 0. 

Several future directions are more or less obvious. We can explore how much the bound 
can be improved by numerically maximizing over differential operators of degree N > 1. We 
can apply the linear functional method to supersymmetric CFT's to get bounds on brane 
tensions in superstring theory. We can try to find linear functional bounds for specific bulk 
CFT's by exploiting knowledge of the bulk spectrum. For example, the most interesting bulk 
universality class for critical quantum circuits is the Monster CFT[15J, which has c = 24 
and Ai = 4. It is interesting because it has no relevant or marginal bulk perturbations. Our 
N = 1 lower bound is g%(24, 4, 1) = 0.0273+. The known conformal boundary conditions |16j 
have g 2 = 1. Preliminary results of numerical calculations for iV up to 37, using the fact 
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that all the bulk scaling dimensions A& are integers > 4, give a bound g 2 > 0.96 which is 
strikingly close to 1. 

The most pressing problem is to overcome the restriction Ai > (c— 1)/12. We expect - 
from consideration of the effective low energy field theory of string theory in the presence of 
branes — that there should be a lower bound on g for all Ai which goes to zero as Aj goes 
to zero. We have shown that in consequence of ^ our present method cannot be extended 
straightforwardly. Some new ideas will be needed. The linear functional method applied 
to the boundary partition function is a practical compromise, well short of the exact lower 
bound that would follow from a complete solution of the conformal bootstrap for boundary 
CFT. We do not know in what direction to improve the linear functional method to get past 
the restriction A x > (c — 1)/12. 
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